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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

AMPLITUDE OF SUPERSONIC DIFFUSER FLOW PULSATIONS 
By William H. Sterbentz and Joseph Davids 


SUMMARY 

A theoretical method for evaluating the stability characteristics 
and the amplitude and the frequency of pulsation of ram- jet engines with- 
out heat addition is presented herein. Experimental verification of 
the theoretical results are included where data were available. 

Theory and experiment show that the pulsation amplitude of a high 
mass-flow-ratio diffuser having no cone surface flow separation increases 
with decreasing mass flow. The theoretical trends for changes in ampli- 
tude, frequency, and mean— pressure recovery with changes in plenum- 
chamber volume were experimentally confirmed. 

For perforated convergent-divergent-type diffusers, a stability 
hysteresis loop was predicted on the pressure-recovery mass -flow-ratio 
curve. At a given mean mass-flow ratio, the higher- value of mean pres- 
sure recovery corresponded to oscillatory flow in the diffuser while 
the lower branch was stable . This hysteresis has been observed 
experimentally . 

The theory indicates that for a ram- jet engine of given diameter, 
the amplitude of pulsation of a supersonic diffuser is increased by 
decreasing the relative size of the plenum chamber with respect to the 
diffuser volume down to a critical value at which oscillations cease. 

In the region of these critical values, the stable mass-flow range of 
the diffuser may be increased either by decreasing the combustion- 
chamber volume or by increasing the length of the diffuser. 


INTRODUCTION 

Supersonic diffusers are usually selected so that at design mass 
flows pulsation-free operation will be obtained. However, there are 
conditions of operation, such as acceleration or maneuvering, which 
may require reduced mass flows in the unstable region. It is therefore 
important that both the conditions under which flow instability will be 
encountered and the resulting pulsation amplitudes and frequency be 
known. 
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Certain resonance criterions which reasonably predicted both the 
conditions for the occurrence of diffuser flow pulsations and the fre- 
quency of pulsation were evaluated assuming that a ram-jet engine acts, 
in effect, as a Helmholtz resonator (reference l) . It was found that a 
diffuser would resonate if the subcritical diffuser characteristic had 
a positive slope of sufficient magnitude . At conditions where the 
slope of this curve was less than the critical value, the ram -jet engine 
would not resonate. A similar criterion for the occurrence of diffuser 

pulsation was suggested from analyses presented in references 2 and 3 . £n 

w 

The occurrence of flow pulsation in external- compression diffusers 
is related in reference 4 to the velocity discontinuity or vortex sheet 
downstream of the intersection of two shock waves . Flow pulsations 
start when this vortex sheet just enters or leaves the diffuser cowl- 
ing. The entrance of the vortex sheet into the diffuser cowling marks 
a change in the slope of the diffuser characteristic which arises from 
the entrance of increasing percentages of low-energy mass flows . The 
amount of change in diffuser characteristic required to induce flow ” 

pulsations in any ram-jet engine may be readily determined from the 
theoretical criterions of reference 1. 

Although certain theoretical criterions as to the nature of dif- 
fuser pulsation have been evaluated, no theoretical evaluation of the 
variables affecting amplitude has been developed. This report, pre- 
pared at the NACA Lewis laboratory, presents a theoretical evaluation 
of the variables which govern pulsation amplitude for the cold-flow 
engine. The theoretical method employed is a refinement of the 
diffuser-pulsation theory presented in reference 1. The method of 
reference 1 assumed small amplitudes of flow disturbances in order 
that quantities such as mass flow and the slope of the diffuser char- 
acteristic could be considered constant. When large finite amplitudes 
are allowed, these quantities must be recomputed in a point-to-point 
calculation throughout a cycle of oscillation. Several numerical com- 
putations which illustrate the nonlinear solutions for pulsation ampli- 
tudes and frequencies as a function of diffuser mass -flow ratio and 
geometry are included herein. Experimental data are presented for 
compaxisqn with the theoretical trends where possible. 


SYMBOLS 

The following symbols and definitions were used in this report: 
A area, sq ft 

a velocity of sound, ft/sec 

f frequency, cycles/sec 
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M Mach number 

m mass flow, slugs/sec 

mo diffuser maximum mass flow, slugs/sec 
Am mass-flow perturbation, slugs/sec 

P total pressure, lb/sq ft abs 

AP pressure perturbation, lb/sq ft 

APj pressure perturbation arising from the inertia properties of the 
perturbed flow, lb/sq ft 

APii pressure perturbation arising from changes in diffuser-shock 
position, lb/sq ft 

t time, sec 

u velocity, ft/sec 

Au velocity perturbation, ft/sec 

V effective storage volume of the ram jet (Vrp - V A ) , cu ft 

V A effective air column volume, cu ft 

Vrp total volume of ram jet, cu ft 

W A mass of the virtual oscillating air column, slugs 

W c mass of gas occupying the effective storage volume of the ram 

jet, slugs 

T ratio of specific heats 

p density, slugs/cu ft 

Subscripts : 
m mean 

0 equivalent free-stream condition 

1 diffuser inlet 

2 diffuser outlet 

4 exhaust -nozzle outlet 
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THEORETICAL ANALYSIS 


The basic concept employed in the analysis of the amplitude of 
pulsation assumes that the resonant properties of a ram-jet engine are 
controlled by the balance between the mass-flow-delivery characteris- 
tics of the diffuser, the storage capacity of the combustion or plenum 
chamber, and the discharge characteristics of the exhaust nozzle. For 
the engine under steady-flow conditions, the mass flow delivered to the 
combustion chamber is equal to the mass flow discharged from the cham- 
ber. The presence of any transient disturbance induces an unbalance 
in the diffuser delivery and exhaust -nozzle discharge rates leading to 
a periodic storage and evacuation of air mass in the combustion chamber 
and concomitant cyclic pressure oscillations . 

A linear differential equation of motion is developed, on the 
basis of the Helmholtz resonator concept (reference l), through a sys- 
tematic calculation of the terms involved in the conservation-of-mass 
equation as functions of the velocity perturbation of the diffuser air 
column and mean-flow parameters of the engine. For conditions of 
unstable operation, this linear vibratory system predicts a continuous 
increase of pulsation to infinite amplitude. In actual systems, of 
course, a positive damping force becomes important at sufficiently 
large amplitudes thus limiting the oscillation. Such a positive damp- 
ing force arises from the nonlinear characteristics of the pressure 
changes incurred from shock oscillations during unsteady operation of 
the diffuser . 

From the following conservation-of-mass equation 


Am 2 


Am 4 + 


cAWc 

dt 


( 1 ) 


the differential equation of motion developed in reference 1 is as 
follows : 


W c W A 


m x rP 2 A l dt 2 


Au 


W c 

"dPg 

T p 2 a 2 M 2 

dm 2 


T-l W A ( m\ a 2 
2 w c (v A 1 2 


_d 

dt 


Au + 


1 _ 

a 2 M 2 



r+1 m 

2r p 2 



0 


( 2 ) 


If second-order effects of terms involving M 2 are neglected, the 
mass flow m x may be stated as 
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“x = 


T p 2 A l M 2 

a 2 


Multiplying by this expression for m x and nondimensionalizing dP 2 /dm 2 
transforms equation (2) to 



Small oscillations were assumed in reference 1 so that the coefficients 
of the differential equation could be treated as constant mean values . 
When large finite amplitudes are allowed, the variation in quantities 

d(P 2 /P 0 ) 

such as — - — r-, P 2 , a-p, and m throughout the cycle must be included. 

d( mg/mo) 

However, equation (2a) can still be applied in a point-to-point calcu- 
lation assuming constant coefficients provided that the increments in 
the cycle are small enough. (The coefficients are, of course, the 
instantaneous values and thus are recomputed for each point.) 

To obtain results for specific diffusers requires a choice of a 
"frozen-state" pressure-recovery mass -flow-ratio curve. The "frozen- 
state" pressure recovery is defined as that recovery which could be 
obtained across the diffuser if steady-flow conditions could be imposed. 
This curve may generally be estimated from the specific diffuser geom- 
etry, the shock configuration, and the diffuser friction losses. 


For a high-mass-flow- ratio single-cone diffuser, peak frozen-state 
pressure recovery will occur at critical mass flow. Pressure losses at 
this condition for such a diffuser will be incurred from a shock con- 
figuration consisting of one oblique shock followed by a normal shock 
plus the estimated subsonic diffuser loss. For lower mass flows, the 
frozen-state pressure recovery decreases because decrements in mass flow 
from the critical condition are marked by the entrance of increasing 
percentages of low-energy mass flows into the inlet as the vortex sheet 
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moves progressively closer to the diffuser centerbody. At zero mass 
flow, the diffuser pressure recovery will he approximately the recovery 
computed across a normal shock at free-stream Mach number. For the 
purposes of this report, it is considered sufficiently accurate to 
define the frozen-state diffuser pressure recovery as varying linearly 
with mass-flow ratio from critical to zero mass flow for the high-mass- 
flow-ratio single-cone diffuser (fig. 1(a)). 

For low-mass -flow-ratio single-cone diffusers, peak pressure recov- 
ery will occur at some value of subcritical mass-flow ratio. The range 
of mass -flow ratios from critical to peak pressure recovery is charac- 
terized by approximately constant total-pressure flow. At further 
reduced mass flows, the frozen-state pressure recovery continuously 
decreases as the vortex sheet enters the inlet and progressively advances 
toward the diffuser cone. By the method presented in reference 5, the 
entrance of the vortex sheet into the inlet and the decrement in pres- 
sure recovery with mass-flow ratio is easily computed for any single- 
cone diffuser. The low-mass -flow curve can again be approximated by 
straight line segments as in figure l(b) . 

Certain cone-type diffusers encounter flow separation from the 
cone surface at reduced mass flows . At present it is not possible to 
predict the mass-flow ratio at which this separation will occur or the 
decrement in pressure recovery it will induce . The method of refer- 
ence 6, however, does predict whether or not a cone-type diffuser will 
encounter flow separation. 

Another type of frozen-state diffuser pressure-recovery mass -flow- 
ratio curve easily established is that curve commonly obtained with 
open-nosed, convergent-divergent, perforated diffusers. For this inlet, 
peak pressure recovery will occur at approximately critical mass flow. 

The theoretical frozen-state pressure recovery varies linearly along a 
line with decrements in mass-flow ratio from critical such as to main- 
tain constant Mach number at the diffuser throat until normal-shock 
pressure recovery at free-stream conditions is reached. This portion 
of the curve (fig. 1(c)) represents the shock in the convergent section 
of the diffuser . The pressure recovery then remains constant at normal- 
shock values with further decrements in mass -flow ratio and this portion 
of the curve represents the normal shock ahead of the inlet . 

After the frozen-state diffuser characteristics are determined, the 
shape of the composite curve of motion and the mean diffuser character- 
istics under oscillating conditions are obtained in the following 
manner: In equation (2a) the following substitutions may be made: 

d 

Au = dt Au 
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Also, the following quantities are assumed to he equivalent: 

w c = P 2 V 

W A = P2 V A 



Thus equation (2a) becomes, after simplification. 



Equation (3) describes a system of curves which relates the variation 
of air-column perturbation velocity Au with the air-column perturba- 
tion acceleration Au. This relation permits the drawing of tangential 
directions in the coordinate system for # the coordinates Au and Au. 
For example, for Au = 0, the slope dAu/dAu is infinite in the sub- 
critical flow region and for Au = 0, the slope dAu/dAu is equal to 
the quantity 



The proper values of pressure, temperature, and mass flow occur- 
ring in the coefficients of equation (3) are determined by successive 
approximation at each point. The mass-flow change given by equa- 
tion (7) of reference 1 stated in terms of Au and Au is 


Am 

“0 



(i) 


(AP X * AP n ) 


+ 


TP 2 Ai 

2 

“ 0*2 


Au 


( 4 ) 


where the pressure changes are 



NACA EM 152124 


^ P 2 V A • 

^1 = -- 27 Au 
a 2 A 1 


or stating P 2 in terms of P 2 ^ m , APj, and APjj 


+ AP )Au 



(5) 


and APjj is equal to the change in frozen-state pressure recovery 
from the mean. The change in a 2 arising from adiabatic temperature 
changes is 


a 2 = a 2 
a 2 a 2 


-ft) 


(?) 


(6) 


where a 2 ^ m and P 2 ^ m are the respective mean values of the diffuser- 
outlet acoustic velocity and pressure. 

An example of the use of equations (4), (5), and (6) follows. It 
is desired to establish the critical mass -flow curve on a Au against 
Au chart for some fixed mean value of mass-flow ratio. The origin on 
this chart, of course, represents the mean mass-flow conditions of the 
problem. The quantities immediately known are Am/m q, P 2 m , APjj, 
a 2,m, geometrical constants, and the variable Au. Quantities to 

be determined are APj, P 2 , a 2 , and the appropriate values of the var- 
iable Au . 

By successive approximations using equations (5) and (6), the val- 
ues of APj, P 2 , and a 2 are easily established for any fixed value 
of Au. (it should be noted that pressure P 2 is the sum of the me an 
pressure P 2 ^ m and the pressure changes APj and APjj.) These quan- 
tities, when substituted in equation (4), permit the evaluation of Au 
and the establishment of one point on the Au against Au chart rep- 
resenting critical flow conditions. When this calculation is repeated 
for other values of Au, the complete critical mass-flow curve is estab- 
lished, Other subcritical mass -flow-ratio curves are established in 
this same manner. 
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The equation defining the supercritical mass-flow curve is 



P 2 m 0 A 1 Au - 


r-i 

2 



= 0 


(?) 


This equation is obtained from a substitution of equations (7), (ll), 
and (12) into equation (lO) of reference 1. For any given value of Au, 
the value of the instantaneous pressure P 2 and the acoustic veloc- 
ity a 2 may readily be determined by successive approximations from 
equation (4). Thus, if the values of P 2 and a 2 for a given value 
of Au are known, the proper value of Au may be computed from equa- 
tion (7) . The repetition of this calculation for other values of Au 
establishes the complete supercritical mass -flow- ratio curve. 

An illustrative solution of equation (3) is plotted in figure 2 
which presents a coordinate system of Au as a function of Au. The 
graph is divided into the subcritical and supercritical flow regions 
by the method just described. Starting at Au = - 600 feet per second 
and Au = 0, the curve goes around in a clockwise direction until super- 
critical flow conditions in the diffuser are established, and there- 
after describes a closed figure continuously. Also, starting from 
rest (Au = 0 and Au = 0) , the same closed curve is reached after a 
number of increasing clockwise spiral turns. 


The variation of Au as a function of time t is obtained from 
a second graphical integration and the results of the example are shown 
in figure 3. The abscissa of a point on figure 2 is represented on the 
ordinate of figure 3 and the ordinate of a point on figure 2 is the 
slope of the curve on figure 3 . Starting from the origin, the ampli- 
tude of Au gradually increases with time until the value of the final 
or natural amplitude is reached. For the case where the initial dis- 
turbance is larger than the final value of amplitude, the pulsation 
diminishes with time to the final amplitude. By use of equations (4), 
(5), and (6) and the graphical variation of Au and Au with time, 
the variation of ram-jet plenum- chamber pressure may be readily 
obtained. 

Some of the effects of changes in engine geometry, flow conditions, 
and diffuser characteristics will be evaluated and compared with avail- 
able experimental data in the following discussion. All theoretical 
amplitude and frequency curves presented were evaluated by the foregoing 
method . 
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RESULTS AND DISCUSSION 
Hi gh-Mas s -Flow-Rat io Dif f us er 

Comparisons of the theoretical and experimental pulsation amplitude 
and frequency for a high-mass-flow- ratio diffuser having no cone surface 
flow separation with two plenum- chamber volumes at a free- stream Mach 
number of 2.0 is presented in figure 4. Schematic drawings of the ram- jet 
engine and diffuser are shown in figure 5. One experimental frequency 
data point was used to establish the proper value of required in 

equation (3) and this point is tailed in figure 4(b). The remainder of 
the frequency curves and the resulting pulsation amplitudes for both 
plenum- chamber volumes were then computed using this value for the vol- 
ume V^. The value of the ratio of the virtual oscillating air column 
to the total ram- jet engine volume V^/Vip for the small plenum- chamber 
engine was 0.740 and for the large plenum- chamber engine, 0.323. 

The theoretical variation and magnitude of pulsation amplitude 
with both volumes is in good agreement with experiment over the range 
of mass flows investigated. The amplitude increased steadily with 
decreases in mass-flow ratio. In both cases, experimental frequency 
decreased gradually with increasing mass-flow ratio compared with the 
nearly constant value predicted by theory. The experimental amplitude 
and frequency decreased as the plenum- chamber volume was increased as 
predicted by the theory. 

From the integration of the pressure-time diagram of a flow pulsa- 
tion ( such as shown on fig . 6) , the average increment or decrement in 
the frozen-state value of diffuser pressure recovery which results from 
the pulsation can be computed. The resulting curves are included in 
figure 4(c) . Flow pulsation caused a decrease in pressure recovery 
below that value assumed for the frozen-state pressure recovery to a 
mass-flow ratio of about 0.61 for the small volume engine and to about 
0.48 for the large volume engine, after which an increase was computed. 

A theoretical effect of an increase in diffuser pressure recovery caused 
by a decrease in plenum- chamber volume was also noted experimentally. 
However, the experimental rise in pressure recovery with mass-flow ratio 
in the range 0.46 < mg/niQ <_ 0.62 was not predicted by the theory. The 

discrepancies may be partly due to the crudity of the assumed frozen- 
state pres sure -recovery mass-flow curve. 

No flow pulsations were found by computation to be sustained for 
supercritical flow conditions and all transient disturbances were found 
to be damped out. 

The theoretical pressure-time wave form only approximately predicts 
the experimental pressure-time fluctuations as can be noted from fig- 
ure 6. Some variations from the experimental values are to be expected 
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because of the simplifying assumptions included in the analysis. A 
closer agreement of the theoretical and experimental wave forms at 
equilibrium by a different analysis are presented in reference 7 . 


Open-Nosed, Convergent -Divergent , Perforated Diffuser 

In the foregoing analysis of the properties of a high-mass -flow- 
ratio diffuser, small point-to-point increments in mass flow, pressure, 
and temperature were used to recompute the point-to-point values of 
the coefficients of equation (3) in obtaining the theoretical solutions. 
As an expediency to conserve time and labor in the following solutions, 
only a few large increments of mass flow, pressure, and temperature 
were used. Thus, although it is expected that if smaller increments 
were employed in the calculation, absolute values of the results would 
be somewhat changed, the general trends will still be the same. 

A theoretical solution of the pulsation characteristics of an 
arbitrary, hypothetical, open -nosed, convergent-divergent, perforated 
diffuser was obtained (see fig. 7). Two modes of operation were found 
to be theoretically possible. In the first mode, if the diffuser were 
started in the zero mass -flow condition with no flow oscillations, 
pulsation-free flow with approximately free-stream normal-shock recov- 
ery would persist over the mass -flow- ratio range to 0.50. At a mass- 
flow ratio of 0.50, however, flow pulsations will set in and the dif- 
fuser pressure recovery will increase several percentage points . Flow 
pulsations, once started, were found theoretically to persist over the 
entire subcritical mass-flow range. This stability hysteresis is 
defined, by theory, to arise from a critical pulsation amplitude such 
that flow perturbations less than this critical value are damped out, 
whereas perturbations greater than this critical value are amplified 
as an oscillatory disturbance to a new natural amplitude. For this 
particular hypothetical diffuser, the critical amplitude is approxi- 
mately the flow disturbance necessary to induce instantaneous values 
of the mass -flow ratio greater than 0.50. The variation of values of 
the critical and the natural amplitude with mass -flow ratio are shown 
on figure 7. 

The existence of a stability hysteresis and the rise in diffuser 
pressure recovery with the onset of flow pulsation are confirmed by 
experiments reported in reference 8. Reference 8 presents the perform- 
ance of a series of perforated, convergent- diver gent diffusers having 
numerous contraction ratios and perforation distributions . Typical 
of the results of this investigation are the data plotted in figure 8 
for the 1.63-e-contraction-ratio diffuser obtained from figure 17 (u) 
of reference 8 . The lower branch of the hysteresis curve was generally 
obtained when the tunnel flow was first initiated. When the mass -flow 
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ratio was increased into the unstable region (from about 0.48 to 0.82), 
the mean-pressure recovery jumped to its oscillatory value. The dif- 
fuser then continued to operate on the upper or oscillatory curve even 
to nearly zero mass flow through this diffuser. 


Some Additional Theoretical Trends 

The effects on amplitude and frequency of varying plenum- chamber 
volume V or diffuser-air-column volume are shown in figures 9 

and 10. In figure 9, the plenum- chamber volume is varied while the 
diffuser volume is held constant. As the plenum- chamber volume is pro- 
gressively decreased, the pulsation amplitude and frequency gradually 
increase until a critical value is reached at which pulsations cease. 

In this region, the stable mass -flow range of a supersonic diffuser may 
be improved by further reductions in plenum- chamber volume. 

In figure 10, the air-column or diffuser volume is varied while 
the plenum-chamber volume is held constant. In this case, as the dif- 
fuser volume is increased, the pulsation amplitude increases and fre- 
quency decreases . Again, a critical ratio of plenum-chamber volume to 
diffuser volume exists such that further increases in diffuser volume 
progressively increase the stable mass-flow-ratio range of the diffuser. 

Thus , from figures 9 and 10, it can be noted that for any ram-jet 
engine of given diameter, the amplitude of pulsation of a supersonic 
diffuser may be decreased by increasing the relative size of the plenum- 
chamber volume with respect to the air-column or diffuser volume. Also, 
the frequency of pulsation for any ram -jet engine of given diameter will 
be decreased as either diffuser volume or plenum- chamber volume or 
both are increased. 


SUMMARY OF RESULTS 

The following summary presents the theoretical and experimental 
results of the effects of engine geometry and diffuser-flow conditions 
on supersonic diffuser pulsation amplitude and frequency. 

1. Theory and experiment show that pulsation amplitude of a high- 
mass-flow- ratio diffuser having no cone surface flow separation ‘increases 
with decreasing mass-flow ratio. The theoretical trends for changes in 
amplitude, frequency, and mean- pressure recovery with changes in 
combustion- chamber volume were experimentally confirmed. 

2. For perforated convergent -divergent -type diffusers a stability 
hysteresis loop was predicted on the pressure-recovery mass -flow-ratio 
curve. At a given mean mass-flow ratio, the higher value of mean-pressure 
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recovery corresponded to oscillatory flow in the diffuser while the 
lower tranch was stable. This hysteresis has been observed experimentally. 

3. The theory indicates that for a ram-jet engine of given diam- 
eter, the amplitude of pulsation of a supersonic diffuser is increased 
by decreasing the relative size of the plenum chamber with respect to 
the diffuser volume to a critical value at which oscillations cease. 

In the region of these critical values, the stable mass-flow range of 
the diffuser may be increased either by decreasing the combustion- 
chamber volume or by increasing the length of the diffuser. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio 
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Mass-flow ratio, mg/mo 

(c) Open- nosed, convergent- divergent, perforated diffuser 


Figure 1 


- Typical variation of frozen- state pressure recovery 
with mass-flow ratio for three diffusers. 
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Figure 2. - Illustrative example of theoretical -growth of minute transient 

disturbance and diminution of large transient disturbance to natural amplitude 
of arbitrary ram jet at fixed mean mass-flow ratio and free-stream Mach number. 
First integration. 


Perturbation velocity, Au, ft/sec 
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0.1 .2 .3 .4 .5 .6 

Time, t, sec 

Figure 3. - Second integration of growth and diminution of transient 
disturbances presented in figure 2. 
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in establishing diffuser-column volume, V A . 



(c) Diffuser characteristic. 


Figure 4. - Comparison of experimental and theoretical effect of 
plenum-chamber volume and mass-flow ratio on pulsation amplitude, 
frequency, and diffuser pressure recovery at free-stream Mach 
number of 2.00. 
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(b) Mass-flow ratio: experiment, 0.41; theory, 0.40. 


Figure 6. - Comparison of theoretical pressure-time variation with 
experiment for high-mass-flow-ratio diffuser (small volume engine). 
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(b) Amplitude. '' 

Figure 9. - Theoretical effect of plenum-chamber volume on diffuser stability, 
pulsation amplitude and frequency. Arbitrary, hypothetical ram- let engine 
with constant air-column volume. 
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(a) Frequency. 



0 4 8 12 16 20 24 

Air-column volume, V^, cu ft 

(b) Amplitude. 


Figure 10. - Theoretical effect of air-column volume on diffuser stability and pulsation 
amplitude and frequency. Arbitrary, hypothetical ram- jet engine with constant plenum- 
chamber volume. 
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